Main sequence stars with asymmetric dark matter 
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We study the effects of feebly or non-annihilating weakly interacting Dark Matter (DM) particles 
on stars that live in DM environments denser than that of our Sun. We find that the energy transport 
mechanism induced by DM particles can produce unusual conditions in the core of Main Sequence 
stars, with effects which can potentially be used to probe DM properties. We find that solar mass 
stars placed in DM densities of p x > 10 2 GeV/cm 3 are sensitive to Spin-Dependent scattering cross- 
section osd > 10 _37 cm 2 and a DM particle mass as low as m x —5 GeV, accessing a parameter range 
weakly constrained by current direct detection experiments. 

PACS numbers: 95.35. +d, 97.10.Cv 



Introduction — Identifying the Dark Matter (DM) is 
one of the most enthralling problems of modern cosmol- 
ogy and particle physics. Weakly Interacting Massive 
Particles (WIMPs) are among the most popular candi- 
dates and are currently searched for with different strate- 
gies, see e.g Ref. [1] for reviews. The effects of WIMPs 
on stars have been recently investigated as possible ad- 
ditional probes for DM searches: in presence of self- 
annihilations, WIMPs captured inside a star can pro- 
vide an exotic source of energy. Whereas the effects on 
the Sun are negligible because of the small DM density 
in the solar neighborhood, they may be remarkable on 
stars living in dense DM environments such as the Galac- 
tic Center, or the first halos to harbor star formation 
at high redshift [2H5]- In general, indirect searches of 
DM would fail for very feebly annihilating DM candi- 
dates. This is the case, for instance, of many asymmetric 
DM (ADM) models (see e.g.jB]), where DM annihilations 
may not occur in presence of an asymmetry between DM 
particles and antiparticles. Still, these particles accumu- 
lated in a star can scatter off nuclei and transport energy 
within its bosom. This effect results in a modification 
of the density and temperature profiles which can lead 
to detectable changes of the solar neutrino fluxes and 
gravity modes [7HT0] ; more compact objects like neutron 
stars may cumulate enough particles to reach the Chan- 
drasekhar mass and drive the collapse of a central black 
hole, [TTJ . Stellar physics can thus provide a strategy 
to indirectly test this class of models. Interestingly, the 
evolution of Main Sequence, solar mass stars placed in 
environments with higher ADM densities has not been 
studied yet (while this paper was being refereed, a re- 
cent study addressing very low-mass stars in ADM has 
been published, [IS])- Here we address this topic, show- 
ing that the energy transport induced by ADM cumulat- 
ing inside such stars can provoke dramatic effects on the 



stellar structure. 

Whether DM particles are actually intrinsically "asym- 
metric" , or if their self-annihilation rate in stellar envi- 
ronment is low enough ((av) 510 _33 cm 3 /s for the Sun) 
is indistinguishable for its effects in stellar evolution, [9]; 
hereafter our definition of "asymmetric" embraces the 
previous condition. 

ADM in stars — Here we adopt the same formalism as 
in [5], to which we address the reader for details. A thor- 
ough description of the underlying theory can be found 
also in [2] and references therein. DM particles can be 
captured by a star via scattering off nuclei, the evolution 
of the total number of DM particles N x inside the star 
reading: 

N X = C-2AN 2 X -EN X (1) 

where C is the particle capture rate over the star, A is 
the annihilation rate and E the evaporation one. For the 
asymmetric DM candidates we are considering, the anni- 
hilation rate is null and evaporation is negligible for the 
physical conditions on which we focus, equation [l] there- 
fore becomes trivial for ADM. Here we recall the linear 
dependence of the capture rate on the spin-dependent 
scattering cross section usd and the environmental DM 
density p x . Captured particles keep on scattering with 
the stellar gas reaching a thermally relaxed distribution 

2/2 

inside the star, n x (r) = n x $e~ r /r x. For the DM masses 
and cross sections we are going to consider, in a Sun- 
like star the thermal DM radius r x is of the order of 
r x <10 9 cm (~ 0.03 i?©), thus making DM particles con- 
fined within the nuclear energy generation region. DM 
particles weakly interacting with the baryons provide an 
energy transport mechanism, additional to the standard 
ones. In their orbits, DM particles scatter in the in- 
nermost regions (called the "inversion core" hereafter) 
absorbing the heat and releasing it in the immediate sur- 
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FIG. 1: The structure of etrans and e nuc for IMq 
star evolved with the ADM parameters /9 x = 10 3 GeV/cm 3 , 
<7SD = 10 -37 cm 2 , m x =10GeV. The profile is taken when the 
central hydrogen fraction, X c =0.2. Also shown for compari- 
son is e nuc of a star evolved without ADM, at the same age. 
The shaded area -innermost of n- shows where etrans has 
negative sign, thus absorbing energy from the stellar gas. 

roudings, within the nuclear energy production region. 

Stars in high ADM densities — The effects of the DM 
energy transport is enhanced for increasing DM densities 
since the number of particles cumulated inside the star 
grows linearly with p x (see the Appendix for a short dis- 
cussion on this). Hence, we have studied the evolution 
of a solar type star (1M , X=0.72, Y=0.266, Z=0.014) 
for increasing DM densities. In this section we fix the 
DM-proton spin dependent cross section crs_D=10 -37 cm 2 
and DM particle mass m x =10GeV. We will discuss the 
effect of varying these parameters in the following sec- 
tions, whereas we keep the stellar velocity through the 
DM halo v=220 km/s. 

We have modified the publicly available DarkStars 
code [14] . in order to use the Spergel & Press formalism 
(Eq. 4.9 in [15]) rather than the Gould & Raffelt one, 
[IS] . Whereas the latter has been proved to correct the 
overestimate -of order unity- present in the previous for- 
malism, we have noticed that in critical conditions such 
as the ones we meet in our study the implementation of 
this formalism may easily induce a unphysical behavior 
of the solution obtained by numerical codes. We have 
found this numerical artifact to show up both with the 
original DarkStars code and the GENEVA stellar evolu- 
tion code [IT], modified for the inclusion of DM effects, 
as we have described in [9]. The effects we describe in 
the following could therefore be overestimated of a factor 
unity, namely they could show up for DM densities p x a 
factor unity bigger than the actually quoted ones. 

We evolve the star from the Zero Age Main Sequence, 
by adopting an environmental /9 x =10 3 GeV/cm 3 , which 
for a NFW DM density profile corresponds to R ga i ~10pc 
from the Galactic center. As long as the absolute value 
of the DM transport energy term etrans is much smaller 



FIG. 2: Stellar (baryonic) temperature profile for two ADM 
parameter set, as from the plot; the solid black line is a model 
evolved without DM. The profile is taken when the central 
hydrogen fraction, X c =0.2. This corresponds at an age of 
approximately 5.3, 5.4, and 5.6 Gyr for the case with no ADM, 
and increasing ADM density p x , respectively. 



than the nuclear energy generation term e nuc , the global 
properties of the star are not modified. We have verified 
that the central temperature and density profile are al- 
tered with respect to the same star evolved without DM, 
according to what we have seen in [9] : neutrino fluxes and 
seismic g-modes of the star get modified with respect to 
the standard case and we defer a systematic analysis of 
this to later studies. 

Here we focus on a much more dramatic effect, taking 
place when etrans becomes comparable with e nuc . For our 
choice of DM parameters m x and <jsd this takes place 
when the total population of DM particles in the star, 
n>tot ~5xl0 47 # (for this particular choice of parameters 
and stellar mass, this takes place approximately 5 Gyr af- 
ter the Zero Age Main Sequence) . In this conditions DM 
is able to transport out of a very central region the en- 
tire energy generated by nuclear reactions, thus resulting 
in an efficient sink of energy for the stellar region within 
the inversion radius fj ~O.O4i?0, namely the region where 
the isothermal WIMP temperature is lower than the local 
baryonic temperature, T x < T\>. Outside of rj, yet within 
the nuclear energy generation region r nuc ^O.1R,0, DM 
particles deposit the energy absorbed from the innermost, 
hotter baryons into the local medium, colder than T x . 
This energy deposit is not crucial: what matters most is 
the effective diminished efficiency of the nuclear energy 
source at the center. The star is forced to compensate for 
this decrease by increasing the nuclear energy production 
outside the inversion core. In Figure[l]we show the radial 
structure of e tra ns compared to e nuc : it is clearly visible 
how the DM energy transport is the mechanism dictating 
the energetics of the region within the inversion radius r^, 
and the energy sink caused by DM in this region induces 
a temperature drop. An example of the resulting temper- 
ature profile is shown in Figure [2] where both the central 
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drop within fj caused by the ADM energy absorption, 
and a raise in the shell between and r nuc are clearly 
visible. Such new temperature profile is the result of 
the new equilibrium reached by the star, which readjusts 
its structure in order to provide the correct energetics; 
the new nuclear reaction rate e nuc is visible in Figure 
[T] These modifications are quite remarkable, modifying 
the chemical evolution of the star as well as its external 
appearance, see more in the next Section. 

The behavior we have described is reached when DM 

particles ctr6 6I10llgll tll&t €-trans ^nuc 

in the central re- 
gion of the star, within r^. For the value of environmental 
DM density p x =10 3 GeV/cm 3 and u=220 km/s adopted 
in this run, this happens quite early during the evolu- 
tion of the star, whereas for DM densities smaller than 
p x =10 3 GcV/cm 3 (with m x =10GeV), the time needed for 
the star to capture enough particles is longer than the MS 
itself, so the mechanism we have described can not take 
place. Conversely, for increasing DM densities p x (and 
therefore the capture rate C) this arrives at shorter times 
during the Main Sequence. 

We have evolved several stellar models for different val- 
ues of the environmental DM density (keeping constant 
(Jsd & n d «i x ), nn ding that remarkable effects on the stel- 
lar structure are present for different DM densities p x . It 
is worthwhile to remark that for increasing p x , DM par- 
ticles start being an effective sink of energy at decreasing 
times. The cumulation of more particles during the MS 
boosts the etrans increasingly, with even more remark- 
able effects: in Figure [3] we show the impact of increas- 
ing ADM densities on the star's evolution. Whereas for 
p x <10 4 GeV/cm 3 \e trans \ > e nuc for R< r nuc , its integral 
over the whole star is smaller than the total nuclear lu- 
minosity of the star, L trans = J \e trans \dV < L nuc : the 
structure reacts to it by finding a new stable configura- 
tion at increased luminosities, by readjusting the nuclear 
reaction rate between and r nuc . On the other hand, at 
very high ADM density (p x JjlO 5 GeV/cm 3 ) the temper- 
ature is mitigated in the whole nuclear energy generation 
region, rather than only in the inversion core. This is be- 
cause the energy reinjected by the WIMPs at R > is 
now sufficient to inflate that region, thus cooling it and 
reducing the efficiency of e nuc . The envelope then con- 
tracts to extract energy from the gravitational potential 
well . 

Summarizing, when L trans < L nuc : the star compen- 
sate for the deficit of nuclear energy generation in the 
center by increasing the nuclear energy generation out- 
side the inversion core. At L tra ns > L nuc : the redistri- 
bution of the energy by the WIMPS is such that e nuc 
is reduced in the whole nuclear core and the star must 
compensate this deficit by the contraction of its envelope. 

Diagnostic power — Is it possible to use such mech- 
anism as a possible probe for ADM parameter space, 
namely: are stellar observables affected from this mech- 
anism at a testable level? Deviations from the standard 
path in the Hertzsprung-Russel (HR) diagram are the 
ideal observable. 



The modifications of the stellar structure described in 
the previous section, induced by ADM, are clearly visible 
in the HR diagram. In our Figure [3] it is possible to see 
how the model described above detaches from a normal 
path in the temperature-luminosity plane once captured 
DM particles have achieved the right number n to t- The 
new distribution of temperature makes the star more lu- 
minous and bigger, thus moving left and upward in the 
plane for low ADM densities. The resulting position of 
the star in the HR diagram is therefore unusual, and more 
so as the star ages and gets toward the end of the MS. 
In Figure [3] it can be appreciated how by increasing p x 
the actual position of the star in the HR changes, finally 
reverting it to lower luminosities as a consequence of the 
dramatic drop in temperature, as explained in the pre- 
vious Section; still, these stars are kept away from the 
usual track. This could indeed be used to identify a pe- 
culiar generation of stars, and probe ADM parameters 
or the distribution of ADM in our Galaxy. However, we 
caution that these observational tests are challenged by 
the difficulty of observing stars close to the galactic cen- 
ter (or in the very center of Dwarf Galaxies, where DM 
environmental conditions could be similar), which makes 
current uncertainties on the position of these stars in the 
HR diagram quite large. 

However, it is worth to remark that, unlike the case of 
annihilating DM, the effects of ADM may be "portable" 
by the star throughout its journey in space. In the case 
of DM annihilating and providing an additional energy 
source to the star, the effects are tightly related to the en- 
vironmental DM density the star is experiencing within 
the short equilibrium time -O(10 6 yr). Self-annihilation 
depletes the DM cumulated inside the star, and a contin- 
uous provision is required to keep the effects going, [2J. 
ADM does not escape from the star: in principle stars 
may capture ADM in a dense environment and migrate 
somewhere else in the Galaxy, and the effects would still 
be visible. 

The effect of varying DM parameters — For the values 
of (Jsd we are considering, the DM energy transport is 
non local, i.e. the DM free path is larger than r x and DM 
particles can efficiently transport energy between distant 
regions in the stellar core. In this case, etrans is en- 
hanced for larger o~sd, until the DM energy transport 
becomes local and so e trans is decreased, this happening 
for (Jsd > 10 — 34 cm 2 . Capture of DM particles is eas- 
ier for lighter ones, and heavier DM particles tend to be 
confined in smaller regions inside the core, both effects 
making the DM transport inside stars more sensitive to 
small values of m x . 

For to x =10 GeV and (7sD=10~ 38 cm 2 the dramatic ef- 
fects on the stellar structure which we have previously 
discussed start to appear for p x >10 5 GeV/cm 3 . This ac- 
tually demonstrates that smaller SD cross sections can 
be explored, although the price to pay is to shrink the 
region of the Galaxy potentially probing this mechanism. 
On the other hand, for £rg£>=10~ 36 cm 2 the same situa- 
tion is reached at lower DM densities, p x =10 2 GeV/cm 3 . 
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FIG. 3: The Hertzsprung-Russel diagram for IMq, for 
varying DM parameters as marked in the plot, with 
<7SD = 10 -37 cm 2 . The solid black square is the starting point 
of our simulation, black solid circles mark an age of 6Gyr, 
not achieved for runs with p x > 10 3 GeV/cm 3 . Black solid 
triangles mark the fall of central hydrogen fraction X c below 
0.1. This takes place at an age of 6.05, 6.36, and 6.39 Gyr for 
the no ADM case, p x = 10 2 ,10 3 GeV/cm 3 , respectively. No- 
tice that the fluctuations may be numerical artifacts, see Ap- 
pendix for details. 

For smaller DM masses (but still above the evaporation 
mass ^5 GeV) it is possible to probe smaller (Tsd, for 
the same value of DM density. 

Effects on massive stars — In principle, the accumula- 
tion of ADM particles inside massive stars should have 
the same effects as in low-mass ones, the transport effects 
of DM not-depending (explicitly) on the burning mech- 
anism (hydrogen via p-p or CNO). However, one crucial 
fact is to be stressed: transport effects start modifying 
dramatically the structure of a star when e trtms is of the 
same order of e nuc within r,;. Two important things are 
to be noticed: i) the stellar luminosity scales non linearly 
with the stellar mass, oc M 3 ' 5 ; ii) the lifetime becomes 
shorter as the stellar mass increases. The latter, com- 
bined with the fact that the capture rate does not scale 
strongly with the stellar mass (roughly C oc M*), leads 
to the conclusion that higher mass stars are less sensitive 
than low-mass ones to the same DM parameters. We 
have evolved models of increasing stellar mass between 
1 and 10 M for the same ADM parameters, that we 
take cr SjD =lCr 37 cm 2 , p x =10 7 GeV/cm 3 , m x =5GeV. The 
5M Q model gets out of the typical path on the HR at 
an age of approximately 90 Myr (X c ^0.2). Stars with 
masses larger than 6Mq evolve normally, i.e. the impact 
of such ADM densities on their evolution is negligible. 

This shows that stars with masses around the solar 
value -or even smaller, as in [T3J- are better ADM probes 
than bigger ones, and yet this is good news, as approxi- 
mately 60% of the stellar mass in our Galaxy is expected 
to be present in the O.IMq <M* <1M range. In the 



very same environment, low mass stars may be affected 
by the dramatic ADM effects, whereas more massive ones 
are insensitive to it. 

It is also worth discussing the effects produced by this 
mechanism on the first stars to form in the Universe, the 
supposedly massive Population III. These objects should 
be naturally placed in very high DM density environ- 
ments, as they are born at the center of a collapsing 
minihalo at high-redhisft, and this is known to enhance 
DM effects onto such stars, [H[S]. 

We have evolved a 100M Q star with metallicity 
Z=10" 4 , with increasing ADM densities, finding that for 
(7s£i=10 -37 cm 2 the effects illustrated throughout this pa- 
per show up at p x >5x 10 10 GeV/cm 3 . Such densities are 
likely to be achieved in the very central regions of a pri- 
mordial star forming halo [18] : yet, recent simulations of 
first star formation show a fragmentation of the proto- 
stellar core, thus yielding to both smaller stellar masses, 
and lower central DM densities [TS] . Whereas no conclu- 
sions can be drawn at the level of the state of the art, it 
is worth to remark that the two effects may compensate 
each other, and PopIII may be indeed affected by ADM. 

Conclusions — In this paper we have shown for the 
first time that feebly or non-annihilating DM, carrying 
weak interactions with the baryons, can have dramatic 
effects on solar mass, Main Sequence stars placed in DM 
environments denser than that of our Sun. This is due 
to the enhancement of DM-driven transport effects, that 
evacuate the nuclear energy produced in center of the 
core of the star. We have shown that once this happens 
this may provoke dramatic changes in the structure and 
external appearance of the star. We have shown that 
such dramatic effects take place in solar mass stars in 
environments with DM densities p x >10 2 GeV/cm 3 , and 
that they can in principle be used as a probe for DM 
particle masses as low as ra x =5GeV and spin-dependent 
cross sections o~sd <^10~ 37 cm 2 . Such parameter values 
are quantitavely comparable with those that could be 
placed by the failed observations of neutron stars in re- 
gions with this very same (fermionic) ADM density, see 
e.g. Figure 1 of [12], yet for as difficult as the observa- 
tion of MS stars can be at the Galactic Center or ADM 
dense environments, they are less challenging than those 
of cold neutron stars. 
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Appendix 

Capture rate — It is to be noticed that we compute 
the capture rate compute for the value v=220 km/s, and 
the p x values we quote in the paper follow this particular 
choice for v (typical of galactic velocity dispersion around 
the solar radius). However, the capture rate grows for de- 
creasing v, and it may therefore be useful to think the 
capture rate in terms of enhancement with respect to 
the solar capture; environments with v lower than the 
one we adopt can in fact be easily found, e.g. above all: 
Dwarf Galaxies. For p x =10 2 GeV/cm 3 , v—220 km/s, the 
enhancement with respect to the solar capture is ^250, 
full formulae can be found in the Appendix of [5]. Self- 
interaction of DM particles can also provide a "boost" 
to capture, in [20] . Numerical issues — In spite of our 
efforts, we have not been able to fully understand the 
reason of the oscillations seen in Figure [3] whether they 
are a physical effect arising from the "bouncing" of the 
central temperature on the WIMP temperature floor, or 
a numerical artifact. We have however checked that the 
existence of oscillations does not affect our results nor 
does change our conclusions. Beside the theoretical con- 
sistency of the intepretation presented in the paper, one 
can be convinced of the actual physical consistency of our 
results with observations based on several numerical ex- 
periments we have performed: a), the deviation from the 
usual HR for mild DM densities takes place before oscil- 
lations start; b), in spite of the variation with timestep of 
the oscillations amplitude (shorter timestep, smaller os- 
cillations), the tracks at varying timestep all converge in 
the same region of the HR; c), no mass loss is associated 
with oscillation, therefore excluding the hypothesis that 
they may drive the evolution of the object by (artificially, 
if they are a numerical artifact) changing its mass. 



